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We present a quantifier of non-classical correlations for bipartite, multi-mode Gaussian states. It 
is derived from the Discriminating Strength measure, introduced for hnite dimensional systems in 
A. Farace et ah, New. J. Phys. 16, 073010 (2014). As the latter the new measure exploits the 
Quantum Chernoff Bound to gauge the susceptibility of the composite system with respect to local 
perturbations induced by unitary gates extracted from a suitable set of allowed transformations 
(the latter being identihed by posing some general requirements). Closed expressions are provided 
for the case of two-mode Gaussian states obtained by squeezing or by linearly mixing via a beam¬ 
splitter a factorized two-mode thermal state. For these density matrices, we study how non-classical 
correlations are related with the entanglement present in the system and with its total photon 
number. 


I. INTRODUCTION 

The presence of correlations in bipartite quantum sys¬ 
tems, that cannot be explained by means of classical 
probability distributions, is one of the main features of 
quantum mechanics. The most important among such 
correlations is surely entanglement [1]. However, re¬ 
cently much attention has been devoted to the study 
and the characterization of correlations that go beyond 
the paradigm of entanglement, being necessary but not 
sufficient for its presence. They are generically known 
under the name of discord-like correlations [2], from the 
name of the first quantifier introduced to measure them, 
i.e. Quantum Discord [3, 4], and appear when the state 
of a (say) bipartite system AB cannot be written as a 
mixture of separable terms that could be locally distin¬ 
guished in one (or both) subsystems. Depending on this 
last choice the set of states with zero Discord is com¬ 
posed by Classical-Quantum (CQ) states, when local dis- 
tinguishability is required for example only on subsystem 
A, or Classical-Classical (CC) states, when it is required 
on both subsystems: 

( 1 ) 

i 

PaP = I*)a (*l ® Ij)i3 01 > (2) 

i,3 

where orthonormal local bases 

(note that differently from the totality of separable 
states, these sets are not convex). 

Since the introduction of the original Quantum Dis¬ 
cord, many other measures M have been introduced to 
quantify the amount of non-classical correlations present 
in a state, approaching the problem from different per¬ 
spectives. All of them, however, must satisfy a set of 
conditions on which the community agrees: 

1. A4 > 0 with equality on CQ (or CC) states; 

2. A4 is invariant under local unitary operations; 


3. A4 reduces to an entanglement monotone on pure 
states; 

with an additional requirement for an asymmetric mea¬ 
sure testing the classicality of subsystem A (i.e. defined 
with respect to CQ states): 

4. A4 is monotonically decreasing under CPTP maps 
on the untested subsystem B. 

A first class of measures uses a geometric description, 
based on the minimization of the distance between the 
considered state and the set of CQ (or CC) states [5-12]. 
Another possible approach is instead rooted on the key 
observation that any state of the form (1) is left invariant 
by a (non-trivial) opportunely chosen local measurement 
or local unitary operation acting on the first subsystem. 
On the other side, if there exists such local action that 
leaves the state unchanged, then the state must be in 
the CQ set. This suggests a way to quantify the asym¬ 
metric non-classicality of a quantum state based on the 
minimum change induced on the whole state by a local 
action, in the form of either a measurement [13-15] or a 
unitary operation. Focusing on the latter, this minimum 
change can be quantified by 

M^a\pab) = min D (pab, Ua) , (3) 

Ua^c> 

where D is some well-behaved non-negative functional 
depending on the state and on a local unitary Ua, and 
5 is a suitably chosen set of unitary operations. Let 
us point out that the idea of characterising quantum 
features of a bipartite system through the detection of 
global alterations induced by local unitary operations is 
not new. For example, in [16] it is shown that global al¬ 
terations due to local cyclic operations are a signature of 
a particular form of non-locality which is beneficial in su- 
perdense coding schemes. However, as already observed 
in [17], the latter investigation is only marginally related 
to the study of discord-like quantum correlations, which 
instead is the goal of the present work. In this respect, 
our assumptions on the set of allowed unitary operations 
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and the subsequent minimization process will greatly dif¬ 
fer from what is done in [16]. 

Expression (3) lies at the heart of several measures of 
correlations known in the literature: the Interferomet¬ 
ric Power (IP) [18] if D is proportional to the Quantum 
Fisher Information [19] of the state when the phase to be 
estimated is encoded by means of Ua, the Discriminating 
Strength (DS) [20] if D is related to the Quantum Cher- 
noff Bound [21] involved in the discrimination between 
the state and its evolution under Ua, and the Discord 
of Response (DoR) [22] if D represents a contractive dis¬ 
tance between such states. The main idea underlying this 
approach is to exploit the effect of the global alteration, 
that in presence of correlations is necessarily induced by 
any local unitary map, in order to have an advantage 
in some information processing task, where not all the 
details characterizing the local evolution are known in 
advance. From this perspective all of the aforementioned 
measures have a clear operational meaning, being strictly 
related to the performance of some tasks: phase esti¬ 
mation for the IP, state discrimination for the DS and 
quantum reading for the DoR. It is also worth stress¬ 
ing that, while being all proper measures of discord-like 
correlations, the quantities introduced in [18-22] are in¬ 
equivalent as no proportionality factor can in general be 
found among them (the different choice of the functional 
D being associated with different potential applications 
of the state pab)- 

Let us emphasize that the minimization set S in Eq. (3) 

plays an important role in guaranteeing that Ma kPab) 
is a proper measure of discord-like correlations. It is 
trivial to see that S cannot be taken to be the whole 
unitary group, otherwise choosing Ua = 1a the corre¬ 
sponding quantifier would be zero for all states. It has 
been shown that condition I. (i.e. Ad = 0 exactly on 
the CQ set) is satisfied for all these measures if the set S 
is taken to be the subset of the unitary operations with 
a fixed non-degenerate spectrum [18, 20]. Therefore one 
can introduce a different well-defined quantifier for each 
spectrum choice. Selecting the optimal one among all 
possible choices is far from a trivial task. The symmetric 
spectrum composed by the cLa = dim "Ha roots of unity 
is conjectured to offer the best resolution, but there is 
no formal proof yet. Eventually, let us point out that al¬ 
though it is possible that other choices of S (apart from 
fixing a non-degenerate spectrum and spanning all possi¬ 
ble basis) might lead to a proper quantifier of discord-like 
correlations, none is known so far. 

Although the previous discussion holds for arbitrary 
systems, most of the quantities involved were first explic¬ 
itly computed only for finite-dimensional systems, often 
qubits. However, there are several cases of interest, for in¬ 
stance in the context of optical interferometry, where one 
might be interested in evaluating the correlations present 
in states of continuous-variable systems, especially Gaus¬ 
sian states [23-25]. In order to answer this question, 
some measures of non-classical correlations have already 
been reformulated by restricting the corresponding op¬ 


erations (measurements or quantum maps) so as to pre¬ 
serve the Gaussian character of the considered states [26- 
28]. However, since this involves a minimization over a 
restricted set of operations, such quantities will in general 
be upper bounds for their original counterparts. Only re¬ 
cently the resolution of the Gaussian minimum-entropy 
conjecture [29, 30] allowed to conclude that the Gaussian 
Quantum Discord [26-28] coincides with the original Dis¬ 
cord on a vast class of states [31]. 

This manuscript is developed within the mindset of 
Gaussian quantum states. In Sec. H we first review some 
basic notions about Gaussian states and Gaussian uni¬ 
tary transformations, then we address the problem of 
characterizing the set S within the Gaussian scenario pre¬ 
viously detailed. We start our discussion with the stan¬ 
dard choice of imposing a fixed spectrum for the unitary 
operations in S with no restrictions on the basis. Next 
we show how the non-degeneracy condition is changed 
when we move from the most general case of arbitrary 
states and operations to the Gaussian scenario. In partic¬ 
ular, if we require some basic properties for the functional 
D {pab, Ua), not all spectra among those associated with 
Gaussian unitary operations can be chosen in order for 
{pab) to be a proper measure of discord-like corre¬ 
lations. Interestingly, the class of allowed sets S for two¬ 
mode systems reduces to a form that has been assumed 
a priori in the specification of some previous measures of 
non-classical correlations to the Gaussian setting [32, 33], 
and that now finds a formal justification. We will then 
focus on the Gaussian version of the DS [20], based on 
the Quantum Ghernoff bound, for which an explicit anal¬ 
ysis in the context of continuous variable systems is still 
missing. In Sec. HI we will obtain a formal expression 
for the Gaussian DS and we will explicitly evaluate it for 
some classes of two-mode states, obtained through linear 
mixing or two-mode squeezing applied on thermal states. 
We will also discuss the relations between Gaussian DS 
and entanglement, squeezing and total number of pho¬ 
tons for these states, drawing our conclusions in Sec. IV. 

II. CHOICE OF LOCAL GAUSSIAN UNITARY 

Let us consider n bosonic modes of a continuous- 
variable system, described by the annihilation operators 
{Qi}r=i- If bfi useful to reorganize them within the 
quadrature vector, satisfying the canonical commutation 
relation: 

f = (ii,pi,... ,x„,p„)T, [f,F]=zr2„, (4) 

where Xi = {hi + dJ)/-\/2, Pi = —i{hi — a\)/\f2. The 
matrix fin in Eq. (4) is the n-mode symplectic form 

= 0 ( -1 0 ) ’ 
i=l ^ ^ 

which enters in the Robertson-Schrddinger uncertainty 
relation that any density matrix p of the system has to 
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satisfy, i.e. 


r -j- i ^ 0, 


of an arbitrary state p there exists a symplectic matrix 
S G Sp{2n, IR) and a set of coefficients such that 


with r being the covariance matrix of p, 


r = Tr[p{f-^,F-^n+]- 


( 7 ) 


r = s 



(14) 


Here {•,•}+ describes the anti-commutator and 

$ = Tr[pf], (8) 

is the associated displacement vector [23-25]. 

In what follows we will focus on the set © formed 
by the Gaussian states, i.e. by the density matrices p 
which can be expressed as Gibbs states of Hamiltonians 
at most quadratic in the quadrature operator f. They 
have a Gaussian characteristic function and are hence 
completely determined by the first and second moments 
of the operator f, i.e. by the vector ^ of Eq. (8) and by 
the matrix F of Eq. (7). Directly associated with the no¬ 
tion of Gaussian states, is the set formed by n-mode 
Gaussian unitary transformations, i.e. by those unitary 
operations which maps 0 into itself. By construction, 
the element of Qn must induce a Heisenberg evolution 
of f which is linear, while ensuring that the commuta¬ 
tion relation (4) holds for the evolved quadrature vector 
operator as well. These two conditions imply that any 
U G Gn can be unambiguously represented as 

U'^rU = lJv+r]u, (9) 

where rju is a vector of and U is a matrix belonging to 
the symplectic group Sp{2n, R), formed by the 2n x 2n 
real matrices S that preserves the symplectic form r2„ 
(5), i.e. 

S e 5'p(2n,R) Srj„ST = (10) 

From Eqs. (7) and (8) one can easily verify that in 
terms of the displacement ^ and the covariance matrix 
r, Eq. (9) induces the following transformations 

+ r^uruT, (11) 

which fully characterize the input-output mapping at the 
level of Gaussian states ©. A special class of Gaussian 
unitaries which will play an important role in the remain¬ 
ing of the paper is provided by the phase-transformations 

[/= ( 12 ) 

with Xj real parameters. In the representation (9) they 
are characterized by a null vector 7}^ =Q and by a block- 
form symplectic matrix 0”^^R(Aj), with R(A) € SO{2) 
being the rotation 

R(A)=f‘=°®^ (13) 

' ^ smA cos A J ' ' 

Let us finally recall that according to the Williamson 
decomposition [23-25, 34] given the covariance matrix F 


The coefficients which are called symplectic 

eigenvalues of F, can be computed as the regular eigen¬ 
values of the matrix |iO„F|, and fulfill the inequality 
Vj > 1 due to constraint (6). In particular this bound 
is saturated on pure states, when all Vj are equal to 1. 
The class of Gaussian states admitting a Williamson 
decomposition with S = 1 are usually denominated 
thermal states. 

We have now all the elements to present our first re¬ 
sult. For this purpose let us assume that the n modes 
of the system are split in two sets: the set A contain¬ 
ing UA modes controlled by Alice and the set B with the 
remaining ub = n — ua modes, controlled by Bob. We 
are interested in constructing a quantifier of the 

form (3), which can be used to characterize the corre¬ 
lations between A and B, by restricting the analysis to 
the case where the allowed initial states pab are Gaus¬ 
sian density matrices of the joint system AB, and un¬ 
der the additional constraint of reducing the set 5 to a 
proper subset of the Gaussian unitary transformations 
Gua^ which operates locally on A. As briefly mentioned 
in the introduction the last assumption is in general not 
justified a priori: analogously to what done in Refs. [26- 
28] it is only motivated by the need of simplifying the 
analysis by forcing also the transformed counterpart of 
PAB under Ua to be Gaussian. We will also assume to 
deal only with sets S composed by all Gaussian unitary 
operations with a certain fixed spectrum, since this is 
the only known choice that in a generic (non-Gaussian) 
framework leads to a measure of discord-like correlations. 
The following Theorem states that under a few basic and 
reasonable assumptions the spectrum characterizing the 
set S must be chosen among those associated with non¬ 
trivial phase-transformations that act locally on every 
mode, i.e. with Gaussian unitary operations as in Eq. 
(12) with Xj not being integer multiplies of 2tt. This can 
be interpreted as the Gaussian counterpart of the non¬ 
degeneracy condition that applies in a generic framework 
where no Gaussian restrictions are imposed. 

Theorem 1. Let us consider a non-negative functional 
D{pabjUa), depending upon a Gaussian state pab o,nd 
a local Gaussian unitary Ua, such that: 

Dl. D{pab,Ua) = 0 PAB = Ua pab U\, 

D2. D{Vb PabVI,Ua) = D{pab,Ua), for all local 
Gaussian unitary operations Vb G Gub , 

D3. D{Va pab GI, Ua) = D{pAB,vi Ua lA), for all lo¬ 
cal Gaussian unitary operations Va G Gua- 


4 


Let us also suppose that given Wa G S, the whole set can 
be written as 

S=\Ua-Ua = VaWaVI K4 e } • (15) 

f S) 

Then the assoeiated quantity \Pab), defined as in 
Eq. (3), satisfies the property 

Ml. M^f^UiocPABUlc) = ■Ma \pab), for all local 
unitary Gaussian operations Uioc = Va®Vb, with 
G Qua Ond B G Qub • 

Furthermore, it will also fulfill the property 

M2. M^^\pab) = 0 PAB = PA® Pb, 

if and only if the operator Wa of Eq. (15) ean he identified 

with a non trivial phase-transformation 
i.e. if and only if the set S writes as 

5{a,} = {Ua : Ua = VAe^^^=^ Fa & Gua} , 

(16) 

where all the real parameters Xj are required not to be 
integer multiples of 271. 

Before going into the details of the proof, let us com¬ 
ment why the assumptions D1~DS, the structural as¬ 
sumption (15), and the properties Ml, M2 appear quite 
natural. 

• The set Dl-Di refers to the “distance-like” quan¬ 
tifier D and are the Gaussian counterparts of anal¬ 
ogous properties which are satisfied by all the 
measures of the form (3) introduced so far, i.e. 
the Interferometric Power [18], the Discriminating 
Strength [20], and the Discord of Response [22]. 
In particular D1 requires that D should be sensi¬ 
tive to all those unitary operations that alter the 

(S) 

considered state. If this were not the case, M 
would not be able to quantify properly the suscep¬ 
tibility of PAB under local unitary maps, on which 
these quantifiers of non-classical correlations rely 
upon. D2 and D3 derive from the reasonable condi¬ 
tion D{pab,Ua) = D{VpABVfiVU aV^), that de¬ 
scribes the independence of the functional D upon 
a change of basis. 

• The hypothesis (15) on the structure of the ele¬ 
ments of S can be interpreted similarly to the re¬ 
quirement of having a fixed spectrum in the finite¬ 
dimensional case [18, 20, 22], without restrictions 
on the basis. 

• Ml, M2 refer instead to the properties that a good 

measure of non-classical correlations is ex¬ 

pected to satisfy. In particular. Ml is just the in¬ 
variance under local unitary operations (see 2. in 
Sec. I) applied to the Gaussian setting, while M2 
is nothing but condition 1. of Sec. I. Indeed, in [35] 


and [36] is shown that the only CQ [see Eq. (1)] 
Gaussian states are those that are completely un¬ 
correlated, respectively in the two-mode and multi- 
mode case. A rederivation of this same fact can be 
found for completeness in Appendix B. The mea- 
sures can therefore be rigorously interpreted 

as Gaussian multi-mode quantifiers of discord-like 
correlations, since they nullify exactly on the same 
set where Quantum Discord does [3, 4]. 

Proof of Theorem 1. The property Ml follows trivially 
from D2 and D3 and from Eq. (15). Let us next 
prove that the form of is necessary to get M2. 

To do so, let us consider an uncorrelated initial Gaus¬ 
sian state Pa® PBi parametrized by a covariance matrix 
Tab = rAOEs and a displacement vector^. Williamson 
decomposition (14), together with Ml, allows us to as¬ 
sume without loss of generality 

UA 

r^ = 0i/,i2, c = o, (17) 

i=i 

with i/j > 1 Vj = l,...,n^. For every product state, 
requirements M2 and D1 impose the existence of a par¬ 
ticular Wa G S that preserves it. From the relation 

wiflT^ = (Wa©1b) f + (18) 

such invariance condition imposes that = 0^^^ and 

( UA \ / riA 

0^.12! Wi= (0 

which, by setting Vj = 1 Vj, in particular implies 

G Sp{2nA) n 0(2n^). (20) 

In Appendix C it is shown that any such Wa can be 
transformed into a direct sum of single-mode rotations 
by means of orthogonal and symplectic matrices, which 
can be adsorbed into the action of Va (appearing in the 
assumed structure of S (15)). In other words we can take 

UA 

Wa=0Ra(A,), (21) 

i=i 

with RA(Aj) G SO{2) identifying hence Wa with the 

Gaussian unitary phase-transformation - 

see discussion above Eq. (13). Finally, in order to have 

a non-trivial measure (condition M2) Wa must be dif¬ 
ferent from the identity in each of the blocks, imposing 
hence all the Xj’s not to be integer multiples of 27r. 

The converse, i.e. proving that property M2 holds 
when assuming in (16), follows from the fact 

that the non trivial phase-transformation e* 

admits the symplectic form (21) and a null displacement 
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vector [see Eq. (13)]. Accordingly, via the correspon¬ 
dence (9) the elements Ua of the set (16) are character- 
ized by symplectic matrices Uyi and vectors r]\j of the 
form 

Ua = ( 0”^! Ra(A,)) (22) 

= (23) 

(A) 

where Va and r]y ' are, respectively, the symplectic ma¬ 
trix and the vector of the Gaussian unitary Va- Now 

let pAB be a Gaussian state which nullifies the quantity 
(S) 

■AdA \Pab)- From D\ this is only possible if there ex¬ 
ists Ua in which leaves such state invariant. At 

the level of covariance matrices this formally implies the 
condition 


Fab = 


Fa 

Foff ^ 

/ UaFaUT 

UaFoff 

^OFF 

Fb y 

\ ^OFF^a 

Fb 


(24) 

where Fab is the covariance matrix of pab expressed in 
blocks form (Fa and F^ being the covariance matrices 
of the reduced density operators associated with the sub¬ 
systems A and B respectively), and where the last term 
refers to the joint covariance matrix of pab after the ac¬ 
tion of Ua- By focusing on the off diagonal-blocks we 
observe that Eq. (24) requires 


V^^Fofb = 0"diRA(A,) V^^Foff, (25) 


for some symplectic matrix Va. Notice that satisfy¬ 
ing the above expression is equivalent to finding, for 
each block, a collection of vectors v solving the eigen¬ 
value equation RA(Aj)'y = v. However, by construction 
Ra(A,) admits as eigenvalues the phases ^ 1. Ac¬ 
cordingly, the only possible solution is i; = 0, i.e. to 
have Vj^^FoFF = 0 in Eq. (25), i.e. Fqff = 0: i-6. to 
have that the Gaussian density matrix pab is a product 
state. □ 


III. GAUSSIAN DISCRIMINATING STRENGTH 

In the following we will focus on a particular quantifier 
of non-classical correlations, the Discriminating Strength 
(DS) introduced in [20], that can be defined using in Eq. 
(3) the following functional: 


Dds{pab, Ua) = 1 — Q{pab, Ua pab U^a)^ 
where Q acts on a pair of states as 

(26) 

Q{Po,Pi)= min Tr [p*p}-'*] . 

sG[0,l] 

(27) 

This yields the following expression for the DS of a bi¬ 
partite state Pab'- 

US^A^ (p) = 1 - max Q (^p, Ua p U\^ 

(28) 

= 1 — max min Tr 

Ua&S sG[0,1] 

P^UaP^-^uI 

.(29) 


The quantity Q that appears in the above definitions 
is the Quantum Ghernoff Bound [21] that intervenes in 
a state-discrimination scenario. We refer to Appendix 
A for a physical interpretation of such quantity, that 
provides the operational meaning of the Discriminating 
Strength measure. 

We will now discuss the problem of obtaining an ex¬ 
pression for a quantifier analogous to the DS for the 
class of Gaussian states. For this purpose we will rely 
on the results of Theorem 1 by observing that the func¬ 
tional in Eq. (26) fulfills the assumptions D1~D3 when 
evaluated over Gaussian states and Gaussian unitaries 
- this can be trivially verified directly from Eq. (27). 
We will hence restrict the optimization set S in (28) 
as described in Sec. II . In particular we face the task 
of evaluating Q{pab,Ua PabU\) for a Gaussian state 
Pab and a Gaussian operation Ua- To do so we can use 
the Williamson decomposition (14), that allows to obtain 
simple expressions for the exponentiated states appear¬ 
ing in the Quantum Ghernoff Bound (27). Indeed, it can 
be shown that, up to a normalization factor 


Tr[piB]=n;UG.(B,), (30) 

p\g is still a Gaussian state whose covariance matrix 
^'ab t)e obtained from (14) just changing the original 
symplectic eigenvalues {i^j}j, pertaining to pab, to the 
set of functions {As{i'j)}j, As and Gs being explicitly 
evaluated as [37]: 


As{x) 

Gs{x) 


(x 0 1)" 0 (cc - 1)" 
{x 0 I)'* — {x — I)"* ’ 
2^ 

{x 0 l)'^ — {x — 1)'* 


(31) 

(32) 


Noticing also that pab and Ua Pab U\ show the same 
symplectic eigenvalues, and exploiting the characteristic 
function formalism [38], the Quantum Ghernoff Bound 
(27) can then be evaluated as [33]: 


Q{pab, UaPabU\) 


min Qse ^ 
se[o,i] 


with 


Qs = 


n- [a,(b,)0Ai,(f,)] 


)/det[: 


A. = 


(s) 

AB 

(«) 


UaF^'-^ 


'u^] 


ab+UaF^V^U^ 


-1 


-(A) 


(33) 


(34) 

(35) 


and with 

Ua = Ua0 1b, = (»?[/\o(^)), (36) 

being respectively the extensions to B [see Eq. (18)] of 
the local symplectic matrix Ua and of the local vector 
r]jj ^ which define the action of the Gaussian unitary Ua 
via the correspondence (9). 

Enforcing hence the restriction (16) discussed in 
Sec. II, in constructing our Gaussian version of the DS 
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we shall use \Ja and 77 ^/^ of the form given in Eqs. (22) 
and (23), optimizing the resulting expression with respect 
to and which parametrize the elements of the 
set Notice, however, that while the minimization 

with respect to can be expected, we can intuitively 
get rid of the displacement parameter r]y . Indeed, when 
dealing with non-classicality measures one is always in¬ 
terested in the worst-case choice of local Ua, that alters 
the states as little as possible. For this reason the choice 
— 0 , leading to a null displacement on the state, 
is expected to be optimal in all situations of interest. 
Exploiting Eqs. (33) and (35), we are now in the posi¬ 
tion of showing that at least for the Gaussian version 
of DS this is indeed the case. For this purpose, observe 
that the matrix -|- is positive definite, 

yielding Ag > 0 with equality for = 0 , and hence 
Qse~^‘ < Qs- This last inequality holds also taking the 
minimum over s and the maximum over Ua, so that we 
can write 


in the Euler decomposition (39), being SO{2) abelian. 
Therefore, Va can be effectively expressed as 

VA{0,x)=-RA{0)S^j\x), (41) 

yielding from Eq. (22) the following functional depen¬ 
dence for the symplectic matrix of Ua, 

VA{e,x) = Ra( 0) sW(cr)RA(A)S^'H-a;)RA(-0), 

(42) 

9 e [0, 27r[ and x € ] — 00 , 00 [ being the parameters over 
which the maximization of Eq. (38) has to be taken for 
fixed A 7 ^ 27rn, n G Z. The resulting GDS for the two¬ 
mode case becomes hence 

QVS^^\pab) (43) 

. [As(^'j)-i-Ai_s(^'j)] 

= 1 — max mm — - =-, 


max Q{p,UapU\)< max Q{p,UapU\). (37) 

Since keeping fixed corresponds to consider a smaller 
maximization set, we can conclude that we are allowed 
to drop all displacements, maximizing only over Va G 
Sp{2nA, R). Therefore, without loss of generality, we fi¬ 
nally define the Gaussian Discriminating Strength (GDS) 
as: 

GVS^a^^Pab) (38) 

= 1 — max mm — =, 

with Ua linked to Va as detailed previously. 


A. Two-mode case 

We will now obtain closed expressions for the GDS 
when referred to special classes of two-mode bipartite 
systems (i.e. n = 2 and ua = bb = 1). In performing the 
optimization of Eq. (38) we use the Euler decomposition 
[39] of a symplectic matrix to parametrize Va, according 
to which every single mode S G Sp{2) can be written as 

S = R( 6 »)S(i)(a;)R( 6 »'), (39) 

where R(0),R(0') G SO{2) as in Eq. (13), and 

a)' <*> 

with X being a real parameter. 

Moreover, since Va intervenes in Eq. (22) always in 
the product Ua = VaRa(A)V)(^, when riA = 1 only 
the squeezing and the rotation on the left are relevant 


where we introduced the function 

F^^He,x) (44) 

= det[r^"]j -f (Ua(0,x) © lB)r^'-*^(UT (0,x) © 1 b)], 

which bares the GDS dependence upon A via Eq. (42). 

We also observe that from the invariance under local 
Gaussian operations of our functional [see property All 
of Theorem 1], the covariance matrix Tab of any two¬ 
mode input Gaussian state pab can be considered in the 
standard form [23, 40, 41]: 

rA = al 2 , rB= 5 l 2 , Fobf = Diag(c, d), (45) 

in terms of 2 x 2 blocks defined as in Eq. (24). Unfor¬ 
tunately even with all these simplifications the explicit 
evaluation of (43) on every two-mode Gaussian state is 
still not trivial. This is because the calculations cannot 
be carried out directly using the coefficients a, b, c, d of 
Eq. (45), since the quantity Qs of Eq. (34) is expressed 
in terms of the Williamson decomposition (14), where an 
high number of parameters intervene in the parametriza- 
tion of the symplectic matrix S [42]. For this reason 
from now on we will consider only two classes of two¬ 
mode Gaussian states, characterized by |c| = |d| in the 
standard form. Such states are obtained from thermal 
states by means of linear mixing (c = d), or two-mode 
squeezing (c = —d). Most importantly, in both cases the 
symplectic matrix S that appears in the Williamson de¬ 
composition (14) of the covariance matrix Tab can be 
described in terms of a single parameter [23, 24]. In par¬ 
ticular, in terms of the Pauli matrix <73 = Diag(l,-1), for 
a two-mode squeezed thermal state Pab identifies 
the symplectic matrix of (14) with 


Ssq(r) = 


coshr I 2 

sinh r <73 

sinh r <73 

coshr I 2 


(46) 
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the connection with Eq. (45) being provided by the ex¬ 
pressions [24, 43] 

vi- V 2 = a-b, {vi + V 2 f = {a + by - , (47) 


where S and C are positive quantities defined as 

S = sinh^(2r), (55) 

C = A\ cosh^(2r) — A^_, (56) 


= {a + lr-A,^ ’ 

vi and V 2 being the symplectic eigenvalues of Tab- For 
a thermal state after a linear mixing via a beam splitter 
instead, Eq. (46) gets replaced by 


Si„,(<^) 


cos (p I 2 

— sin p I 2 

sin p I 2 

cos p I 2 


(49) 


for the squeezed thermal states of Eq. (46), and 

S = A^_siT?{2(j)), (57) 

C = A\- A^_cos^{2(j)), (58) 

for the thermal states after a linear mixing of 

Eq. (49). Notice that in Eq. (54) there is no dependence 
upon 9 and that the left hand side reaches the minimum 
for X = 0. Accordingly we can write 


while Eqs. (47) and (48) by the identities 

vi + V 2 = a + b, {vi — 1 / 2 ^ = {a — b^ -|- 4c^, 


sin^(2(/)) 


4c^ 

(a — by -I- 4c^ 


(50) 

gvs^A^ 

(pS) 

(51) 

gvs^y'’ ( 

>s) 


sinh^(2r) sin^(A/2) 



+ 1 1 

5^ ^ 

5 

-I- sinh"(2r) sin" (A/2) 
(59) 

2p) sin" (A/2) 


r-‘l 

-1- sin" (20) sin"(A/2) 


1. Explicit evaluation on linear mixing and two-mode 
squeezing of thermal states 


It is important to stress that the minimum over the 
parameter s in Eq. (38) is in general difhcult to evalu¬ 
ate, because it has to be performed for a generic choice 
of the parameters 9 and x which, via Eq. (42), define 
the symplectic matrices over which we have to take 
the minimization. A lower bound on GDS can be ob¬ 
tained setting s = 1/2 throughout the computation. In 
[33] it is shown that when the rotation parameter in (42) 
is set to be A = ±7r/2, this bound is achieved for every 
Gaussian state. In Appendix D we show that for the two 
aforementioned classes of states, such minimum is actu¬ 
ally reached in s = 1/2 for every choice of A. However, 
we have numerical evidences that this is no more true 
for a generic two-mode Gaussian state not fulfilling the 
symmetric condition jcj = jdj. 

Applying the aforementioned result in Appendix D, we 
now evaluate (43) setting directly s = 1/2, obtaining 


gvsfipAB) = 1 - 


MAl-A^) 


(52) 


where in writing the nominator we introduced the quan¬ 
tities 


Ai/2(*^i) ± Ai/2{v2) 


(53) 


After long but straightforward calculations, we can also 
express the term at the denominator as 

{9, x) = [4(A^ - A^) + 4 sin"(A/2) S ]" 

-I- 16sinh"(2x) sin"(A)(A^ — A" ) C, (54) 


(60) 

The optimal measure for these classes is therefore the one 
described by the parameter A = tt, that yields the max¬ 
imum amount of correlations for every given state (e.g. 
see Eigure 1). Notice also that due to the dependence 
upon A^, the obtained result is invariant under the ex¬ 
change of the two subsystems [see Eq. (53)]. This feature 
appeared also in the Gaussian Interferometric Power [32] , 
and it is a peculiarity of the considered classes of states. 
Qualitatively we can see that fixing the squeezing r or 
the linear mixing parameter (p, GDS approaches its max¬ 
imum value 1 when A?. ~ A^, that from (53) corresponds 
of having a big gap between the symplectic eigenvalues vi 
and V 2 of Tab (i-O. in the excitation numbers of the ther¬ 
mal states). This is analogous to the result pointed out 
in [33] for the Gaussian Discord of Response, where the 
amount of correlations can be increased with the asym¬ 
metry on the number of thermal photons between the two 
subsystems. Notice also that symmetric states {vi = 1 ^ 2 ) 
yield the same amount of non-classicality as pure states 
(which is zero for the linear mixing case), and that for 
two-mode squeezed states the maximum value of 1 can be 
obtained for every thermal state in the limit of r —> 00 . 
Similarly, the maximum amount of correlations obtained 
through linear mixing is achieved by means of a balanced 
beam splitter {/p = 7r/4), but it can approach 1 only for 
highly asymmetric initial thermal states. For more infor¬ 
mations on the correlating power of beam splitters, see 
[44-46]. 

The expressions obtained here are written in terms of 
the parameters that intervene in the Williamson decom¬ 
position (14): exploiting Eqs. (47), (48), (50) and (51) we 
can, however, convert them in terms of the the parame¬ 
ters a, 6, c and d which appears in the standard form (45). 















GDSa*''’ 



FIG. 1. (Color online) GDS of the two-mode squeezed vacnnm 
state as a function of the squeezing parameter. In the plot 
we have used a set of exponentially decreasing parameters 
A = ^. From left to right we find the curves for fc = 0,1,... 6, 
the optimal choice being associated with k — 0, i.e. A = tt. 


In particular for symmetric two-mode squeezed thermal 
states, i.e. a = b or vi = V 2 , this yields 




sin^(A/2) 
cos^(A/2) ’ 


(61) 


that when A = -n 12 coincides with the related re¬ 
sult on the Quantum Chernoff Bound obtained in [33]. 
For asymmetric states, however, the contribution of 
(A_/a 1+)^ becomes cumbersome when expressed as a 
function of standard form parameters. An expression 
analogous to (61) can be derived also for the linear mix¬ 
ing case. Indeed for a = b (which in this case does not 
correspond to have vi = 1 ^ 2 ) the GDS reads: 




(62) 


4c^ sin^(A/2) 


{a + c)^ — 1 + y/{a — — 1 — cos2(A/2) 


B. Relation between GDS and entanglement - total 
number of photons 

In the previous sub-section we showed that keeping the 
squeezing or the mixing parameter fixed, the amount of 
correlations increases with the asymmetry in the number 
of thermal photons used. Here we want to discuss how 
the GDS behaves for the same two classes of states with 
respect to entanglement or to the total number of pho¬ 
tons. In other words, our goal here is the identification of 
the optimal state achieving the maximum amount of cor¬ 
relations, when we fix either the degree of entanglement 
or the total energy contained in the state. 

We begin by considering entanglement on the squeezed 
class only, being always separable. To quantify it, 
we can use the logarithmic negativity [47, 48] , defined as 

£■ = max{—log(z/_), 0}, (63) 


where v- is the smallest symplectic eigenvalue of the par¬ 
tially transposed state. Exploiting the method based on 
the symplectic invariants [42, 43], can be explicitly 
written as: 


V- 


A - VA2 -ddetr^B 


(64) 


where A = det F^i-l-det Fs—2 det ro_FF- For a two-mode 
squeezed thermal state this corresponds to det F ab = 

v\v2 


A = cosh(4r) 

so that an explicit expression for £(vi,V 2 ,r) can be ob¬ 
tained (see Eq. (E5) in Appendix E if interested). In 
particular, for the pure two-mode squeezed vacuum one 
has £(I,I,r) = 2jrj. Therefore, for this class of pure 
states GDS is a monotonic function of the logarithmic 
negativity. At least in this particular case, we can see 
how property 3. of a good measure of non-classical cor¬ 
relations (see Sec. I), namely of being an entanglement 
monotone on pure states, is satisfied. 

We numerically evaluated GDS (with the optimal value 
A = tt) and the logarithmic negativity for a set of 10® two¬ 
mode squeezed thermal states randomly generated with 
the constraints 1 < = 1/2 < 20 and 0 < r < 5. The re¬ 

sult is shown in Figure 2, where we put in evidence the be¬ 
havior of symmetric states with vi = V 2 in both the pure 
(red lower line) and the mixed case (black dashed lines). 
Notice that for every fixed value of entanglement GDS is 
minimized on the pure two-mode squeezed vacuum, on 
which the entanglement is easily related to the squeezing 
parameter by the relation £(I,I,r) = 2jrj. Such bound 
holds even if we allow vi to be different from 1/2 ■ On the 
other side, considering mixed states the value of GDS is 
allowed to go as close to I as desired. These facts are for¬ 
mally stated in the following proposition, which is proven 
in Appendix E. 

Proposition 1. Consider Q'DS'^^\vi,V 2 ,r) for a two¬ 
mode squeezed thermal state, and its logarithmic negativ¬ 
ity £{ui, 1 ^ 2 , r). Then: 

• Fixed £, the GDS is minimized on pure states: 

gvs^^\iz,,,z 2 ,r) > {l,l,£{v^,iz 2 ,r)/ 2 )- 

• For all £ > 0 and A such that 

gVS^^\l,l,£/2.)<^<l, (66) 

there exists a symmetric state with izi = 1/2 = iz 
with entanglement £ and gDS^^^ = A. 

This means that if entanglement is the quantity in 
which we are interested in optimizing, then pure states 
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FIG. 2. (Color online) GDS with respect to logarithmic 
negativity £ for a set of 10® randomly generated two-mode 
squeezed thermal states with 1 < = V 2 = v < 20 and 

0 < r < 5. The curves corresponding to pure states (red, 
lower bound) and to symmetric states with iz = 2,3,4, 5, 6, 7 
(black dashed lines from bottom to top) are also shown. The 
plot is realized with A = tt. 

offer the worst performance, whereas mixed states with 
the same entanglement yield always a greater GDS. 

The situation, however, changes if we take into account 
also the energy that must be used to produce such ther¬ 
mal excitations. This can be done comparing GDS with 
the total number of photons N of the state. For a linear 
mixed thermal state this reads: 


FIG. 3. (Color online) Gaussian DS with respect to the total 
number of photons for a set of 10® randomly generated two¬ 
mode squeezed thermal states with l<u\=V 2 = v<S: and 
0 < r < 7. The curves corresponding to pure states (red, 
upper bound) and to symmetric states with iz = 2, 3,4, 5, 6, 7 
(black dashed lines from left to right) are also shown. The 
blue dot-dashed line that crosses the others corresponds to 
the optimal GDS achieved only by means of linear mixing of 
thermal states. The plot is realized with A = tt. 

In order to do an equivalent calculation for a two-mode 
squeezed state characterized by the symplectic matrix 
(46), we need to consider a thermal state evolved with 
the operator S'sq(r) = that modifies the total 

number of photons as: 


Mm 


Tr 



y^h)S\^{4>)PthSl^{4>) 


Sl^{a^a + Ph)Ssci = 2sinh^ r-|-cosh(2r)(aM-|-S^6). (70) 


= Tr 


(oM -I- y^h)p 


'th 


Vi + V2 


- 1 , 


(67) 


Therefore in this case the total number of photons de¬ 
pends also upon the squeezing: 


where we used the fact that a beam splitter does not 
change the total number of photons. In the previous ex¬ 
pression Pth is file thermal state with symplectic eigenval¬ 
ues vi, h' 2 , while Sini{4>) = is the unitary op¬ 

eration that corresponds to the symplectic matrix Si„j(0) 
of Eq. (49) through relations (9) . In this case the total 
number of photons fixes exactly the sum vi + 1 ^ 2 , there¬ 
fore the maximum of GDS (60) is obtained when the 
gap between them is as big as possible: vi = 2Mm + 1 
and i /2 = 1 or viceversa. From this it follows that the 
maximum amount of correlations that a balanced beam 
splitter can generate, starting from thermal states with 
a total number of photons N, is given by: 




N — N cos A 
2 + N — N cos A 


( 68 ) 


With the optimal choice A = tt, it reduces to the simple 
expression 


GDS‘:^'> 


/ (lm,opt)x 

\Pab ) 


N 

N + V 


(69) 


that is plotted as a blue dot-dashed line in Figure 3. 


Mq = Tr 


(oM -I- Pb)Ssq{r)pthSlq{r) 

= 2 sinh^ r + cosh(2r)Tr (SM + Pb)pth 
= cosh(2r) ^ — 1, 


(71) 


For the case of pure states this expression reduces to 
fVsq(l,l,r) = 2sinh^r. Analogously to what have been 
previously done for entanglement, we can numerically 
study the relation between GDS (with A = tt) and the to¬ 
tal number of photons for a set of 10® randomly generated 
Gaussian states with 1 < bi = 1^2 8 and 0 < r < 7. 

The result is plotted in Figure 3, where the behaviour of 
pure states (red upper line) and symmetric mixed states 
with h'l = V 2 (black dashed lines) is evidenced. Notice 
that, differently than before, pure states corresponds to 
the GDS upper bound for of a two-mode squeezed ther¬ 
mal state, when the number photons is kept fixed. As in 
the entanglement case, such bound holds also for z/i ^ V 2 - 
Every other value of correlations below the threshold of 
pure states can instead be achieved. Both these facts 
are formally stated by the following proposition, that is 
proven in Appendix E. 
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Proposition 2. Consider Q'DS^^\vi,V 2 ,‘>') for a two¬ 
mode squeezed thermal state, and its total number o/p/io- 
tons Nsq(yi,V 2 ,r). Then: 

• Fixed Ngq, is maximized on pure states: 

Q'DS^^\i'i,v2,r) < 

GDS'P (” 1 , 1 , ] ; 

• For all Nsq > 0 and A such that 

0 < A < ^1,1, arcsinh Nsq/2^'^ , (72) 

there exists a symmetric state with = 12 ^ = r' 
with total number of photons Nsq and = A. 

Comparing Proposition 1 with Proposition 2 we can see 
that, depending on which resource we take into account, 
pure states can be either the optimal or the worst possible 
choice. They indeed allow to reach the maximum amount 
of correlations when the total number of photons is fixed. 
However, if thermal excitations are already present in the 
system, maybe due to some noise, then the mixedness 
that is introduced increases the non-classical behaviour 
of the state when the same squeezing parameter or the 
same amount of entanglement is considered. 

Eventually, in Figure 3 one can see also a comparison 
between the maximum CDS that can be activated from 
thermal states, by means of a balanced beam splitter, or 
via a two-mode squeezing operation leading to the same 
total number of photons. It can be appreciated how the 
two-mode squeezed vacuum always performs better, but 
at the cost of utilizing highly-squeezed light, which is 
much harder to produce with respect to the combination 
of thermal light and linear mixing through a balanced 
beam splitter. 

IV. CONCLUSIONS 

The susceptibility of a state to the less disturbing lo¬ 
cal unitary map has been widely used in the literature to 
quantify the amount of non-classical correlations. In or¬ 
der to obtain a bona-hde measure, which nullifies only on 
the set of CQ states, such unitary is usually chosen among 
those operations characterized by a fixed non-degenerate 
spectrum. In the first part of this paper we studied how 
such condition is modified when a Gaussian restriction 
is imposed on both the state and the unitary maps in¬ 
volved in the minimization. In particular, in Sec. II we 
showed that, under few very reasonable assumptions, all 
local unitary operations must be chosen among those that 
can be obtained from a non-trivial operator of the form 
by means of local Gaussian unitary maps. 
Although its single-mode version has been assumed sev¬ 
eral times in the literature, to the best of our knowledge 


no explicit proof of its necessity was known so far. More¬ 
over, we discussed how usually the optimization over the 
remaining local Gaussian unitary map can be performed 
only by considering its action on the covariance matrix, 
without taking into account the possibility of a displace¬ 
ment. 


After this general discussion, we focused on a partic¬ 
ular measure of non-classical correlations: the Discrim¬ 
inating Strength. With a clear operational meaning, it 
aims at enlightening the usefulness of correlations in a 
state discrimination setup. We obtained an expression 
for its Gaussian version (GDS), defined considering Gaus¬ 
sian input states and restricting the optimization set S 
to a subset of all Gaussianity-preserving unitary maps, 
as described in the first part of this manuscript. We then 
computed the obtained quantity on two simple, but still 
relevant, classes of two-mode states: those obtained from 
thermal states by means of a linear mixing or a two-mode 
squeezing operation. Explicit expressions for GDS have 
been provided in these cases, in terms of the parameters 
intervening in their Williamson decomposition, as well 
as in terms of their standard form defined up to local 
unitary maps. We showed that given an initial thermal 
state, the maximum amount of non-classical correlations 
that can be obtained by mixing linearly the two modes 
is achieved with a balanced beam splitter, whereas in¬ 
creasing the squeezing parameter always leads to a big¬ 
ger amount of non-classicality. Furthermore, asymmetric 
states with a different number of thermal excitations in 
the two modes, i.e. with symplectic eigenvalues vi V 2 , 
are always more correlated than their symmetric coun¬ 
terpart. 


We compared GDS in a two-mode squeezed state with 
the amount of entanglement, as measured by the log¬ 
arithmic negativity. We showed that, for every fixed 
value of entanglement, pure states offer the worst per¬ 
formance in terms of amount of correlations, whereas 
mixed states can lead to values of GDS as close to the 
maximum bound of 1 as desired. However, this is due to 
the fact that states with initial strong thermal popula¬ 
tions can be highly squeezed despite their small amount 
of entanglement. Therefore, it is also interesting to com¬ 
pare different states at fixed energy. From the analysis of 
the relation between GDS and photon number emerged 
that when the total energy is kept fixed the pure two¬ 
mode squeezed vacuum is indeed the option leading to the 
biggest amount of correlations. Eventually, we showed 
that mixing different thermal states in a balanced beam 
splitter still allows the generation of non-classical correla¬ 
tions. Despite the fact that this scenario does not achieve 
values as high as a two-mode squeezed vacuum, it has the 
advantage of not requiring highly-squeezed light, which 
is harder to obtain experimentally with respect to many- 
photons thermal states. 
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In this section we review the state-discrimination pro¬ 
tocol that intervenes in the definition of the measure of 
non-classical correlations known as DS introduced in [20] . 
Before going into details, it is worth to recall the Quan¬ 
tum Chernoff Bound (QCB) [21], that describes the prob¬ 
ability of error when discriminating between two known 
quantum states po and pi, available in M copies. Opti¬ 
mizing over all Positive Operator Valued Measurements 
(POVM), if there is a 50% a priori probability of having 
pf^ or pf^ one has [49]: 

P^e^r!rmniPO,Pl) = ^ (l " \\pf^ " pf""l|l) , (Al) 


corresponding to an optimal strategy that discriminates 
between the positive and negative eigenspaces of the dif¬ 
ference — pf^. In the limit M 3> 1 such probability 
scales exponentially in the number of copies: 

(A2) 

where we defined 


^{Po,Pi) 


lim 

M—¥00 


log ^err^min (Po,Pi) 


M 


(A3) 


Such limit intervenes in the definition of the QCB [21], 
that appears in Sec. Ill of the main text: 


i(Po^Pi) = Q(po,pi) = niin Tr [pgpj '*] . (A4) 

sG[0,l] 


Notice that in particular the QCB can be bounded as: 


0 < Q(po,pi) < Tr 




< 1 . 


(A5) 


The DS is based on a quantum illumination scenario 
[50-53], where one is interested in obtaining informa¬ 
tions about an environment by probing it with a bipartite 
state: the actual probe A and a reference i?. In particular 
the DS considers a situation in which one is interested in 
detecting the possible application of a local unitary op¬ 
eration performed by a third party, that is only partially 
cooperative. The protocol can be thought as follows (see 
Figure 4): Alice and Bob choose the initial probing state 
PAB, and produce M copies of it. Alice’s local state is 
then sent to a third party Charlie, who can apply or 
not (with 50% probability) a unitary operation Ua of his 
choice. Charlie then communicates Ua (without saying 
whether he applied it or not) and gives back to Alice the 
states, on which she can perform a generic joint POVM 
on the whole system AB with the help of Bob. Their aim 
is to discriminate between pab and UaPabU\, in order 
to detect whether the local rotation actually took place 
or not. 

Notice that the chosen Ua, not known at the state 
preparation stage, is eventually revealed in order to al¬ 
low Alice and Bob to perform the best possible measure¬ 
ment, whose outcome can be described by the QCB. This 



FIG. 4. Scheme for the state discrimination protocol on which 
the Discriminating Strength is defined. Alice and Bob prepare 
n copies of a state, without knowing which local unitary oper¬ 
ation Ua will be applied in the following. Alice’s subsystems 
are then modified locally by a third party (Charlie), that can 
choose with the same probability to: (i) leave them unchanged 
or {ii) to apply a local rotation Ua G 5 of his choice. The 
picked unitary is then communicated to Alice and Bob, that 
can use this information to perform the optimal joint mea¬ 
surement to detect whether Charlie chose option (i) or option 
{ii). 


guarantees that the worst-case success probability only 
depends on the initial chosen state pab and on the set 
S of the allowed unitary operations Ua- As can be ex¬ 
pected, the presence of non-classical correlations in the 
initial state pab turns out to be useful. Consider for 
example the structure of a CQ (uncorrelated) state (1): 
indeed, if Charlie chooses a unitary operation diagonal 
in the same basis {|*)^}i in which the CQ state is ex¬ 
panded, then the whole state pab is left unmodified and 
Alice and Bob have no way to detect the rotation. The 
more the initial state pab is susceptible to local alter¬ 
ations, the higher the probability of success of Alice and 
Bob will be. Therefore, the Discriminating Strength of 
the initial state pab defined in Eq. (28) of the main text 
quantifies their worst performance in such discrimination 
task. 

In a previous work [20] it was shown how the DS is a 
bona-fide measure of non-classical correlations, satisfying 
the four properties stated in Sec. I, when S is chosen to 
be the set of all local unitary operations with a hxed and 
non-degenerate spectrum. 


Appendix B: CQ Gaussian states are completely 
uncorrelated 

In this appendix we rederive that only completely un¬ 
correlated Gaussian states, i.e. decomposed s\s pA® Pb, 
are CQ (i.e. with a zero Quantum Discord, see Sec. I), 
result that has been hrstly obtained in [36]. The discus¬ 
sion will follow the ideas of [35] , where only the two-mode 
case is considered. 

The core of the proof relies on a result involving 
informationally-complete POVM (IC-POVM) [54, 55]. A 
POVM {Mk}k, with Mfc = 1, is said to be complete 
if it allows to reconstruct the state, i.e. if there exist 
operators such that a = VfcTr[(TMfe], for every 
density matrix a of the system being measured. If we 
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perform an IC-POVM on the subsystem i? of a bipartite 
state PAB^ we are left with the conditional states: 


PA\k = 


TrB[pA5-^fc] 

Pk 


Pk =Ti[pABMk]. (Bl) 


Rahimi-Keshari et al. in [35] have proven the following 
statement. 


Proposition B.l. The necessary and sufficient condi¬ 
tion for a bipartite system pab to be CQ is that, for every 
IC-POVM {Mk\k performed on subsystem B, the condi¬ 
tional states {pA\k}k commute with one another, i.e. 

[PA|fc,PA|fc'] = 0, yk,k'. (B2) 


Based on this result, in order to prove that a CQ state 
p S 0 has to be completely uncorrelated, we can just 
prove that if its covariance matrix has some correlations 
PoFF 7 ^ 0 [see (24)] then there exist an IC-POVM lead¬ 
ing to two conditional states PA\k-i PA\k' that do not com¬ 
mute. In the following we will consider a particular IC- 
POVM for continuous variable systems: the multi-mode 
heterodyne detection, that projects on coherent states 

= (B3) 

with P G C"®. Notice that the completeness follows from 
the correspondence between quantum states and quasi¬ 
probability distributions on phase space, that can be re¬ 
constructed from such measurement. 

To describe the conditional state pA\p after a hetero¬ 
dyne measurement on subsystem B, we will use the fol¬ 
lowing well known result on Gaussian POVMs [24, 56]: 

Proposition B.2. Let us consider a bipartite Gaussian 
state Pab described by a covariance matrix and a dis¬ 
placement vector given by: 


^ AB = 


Pa 

Boff 

Boff 

Bb 



whose subsystem B is projected on a Gaussian state de¬ 
scribed by Fjji and After the measurement, the state 
left in A has covariance matrix and displacement vector 
given by: 

r^lm = —- Pq^^, (B5) 

i i m 

- r OFF ^ p (Cb - ■ (B6) 

i- B + i-m 

In our case, up to local Gaussian operators (that do not 
alter the amount of correlations), we can always choose 
^A — ^B — 0- Moreover, in the heterodyne detection 
Pm = 1 b and 0 G The conditional state pA\p 
will be therefore characterized by 

r.4|P = Ba — PoBF:p; —j ^OFF^ (^7) 

^A\P = ^OFF^ - 

1 B + 1b 


Notice that being Pa |/3 independent on /3, all the con¬ 
ditional states pa\p differ only for their first moments. 
Since Poff 7 ^ 0 by hypothesis and being P^ -f 1^ G 
GL{2nB), there exists a vector /8 q such that 7 ^ 0, 

whereas choosing /3 = 0 one has = 0. In the follow¬ 
ing we will focus on the couple of states pa\o and PA\Pg, 
proving that they do not commute. 

To do so, we will use the fact that the product pip 2 of 
two states pi,p 2 & ® still has a Gaussian characteristic 
function (possibly with non-real covariance matrix P 12 ) 
[38]: 


Xp,p.(w)=x(0)e-J“^^-“-*«-“. (B9) 

In particular, can be written as [38]: 

ei2=^l-(ri-*ll):p4T^(^l-e2), (BIO) 

1 1 -h 1 2 

where and Pi, i = 1,2, are respectively the displace¬ 
ments and covariance matrices of pi, p 2 and Ft is the sym- 
plectic form (5). Condition [pi, P 2 ] = 0 implies 
because to different characteristic functions correspond 
different quantum states. From (BIO) and its counter¬ 
part for ^ 2 ii such equality is equivalent to = ^ 2 - 
In conclusion, given pab S © with Poff 7 ^ 0, we 
found a continuous variable IC-POVM (heterodyne de¬ 
tection) yielding two conditional states PA\Pg and pa\o 
characterized by different displacements Ca\Pq 7^ ^A| 0 ' 
Since this implies their non-commutativity, from Propo¬ 
sition (B.l) it follows that pab cannot be CQ. Conversely, 
if a Gaussian state pab = Pa® Pb is completely uncor¬ 
related, it can be explicitly written in a CQ form (1) by 
diagonalizing its local components pA and ps- 


Appendix C: Standard form for a symplectic 
orthogonal matrix 

It is a well known fact of linear algebra that a special 
orthogonal matrix W G SO(2n) can be reduced to a 
direct sum of n two-dimensional rotations (40) by means 
of an orthogonal change of basis: 


W = M 


0R(A,) MT, 




M G 0(2n). (Cl) 


In this appendix we will prove that if W is also symplec¬ 
tic, then the same decomposition can be obtained with 
M G Sp{2n, R) n 0{2n). From now on it is convenient to 
adopt a different quadrature operator vector f, in which 
the position-like operators come before the momentum¬ 
like ones. With this definition all the relevant properties 
hold with respect to the symplectic form 


J = 



(C2) 
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with zeros in the diagonal blocks. Notice that with this 
different ordering the standard form of W, equivalent to 
(Cl), will be 


= 


0, cos Xj 

0,-sinAj 

BjSinXj 

0j cos Xj 


(C3) 


In order to complete our proof, we will need two general 
properties of the symplectic group [57], that are stated 
in the following lemma. The proof is sketched for com¬ 
pleteness. Note that a complex matrix Me is said to be 
symplectic if it satisfies the same formal requirement as 
in the real case: 


MeJMJ = J, 


and we will write Me S Sp{2n, C). 
Lemma C.l. Let us consider the matrix 


1 

' 1„ 

i\n 

V2 


In 


(C4) 


(C5) 


that given a real matrix M allows for the definition of its 
complex correspondent as 

Me = KMK-\ (C6) 


Then: 


(i) M G S'p(2n,R) ^ G Sp{2n,C); 

(ii) M G Sp{2n, R) n 0(2n) ^ = U 0 U% 

for some U G U(n). 

Proof (i) The proof of the first statement follows trivially 
from the symmetry of K and the observation 

KJK = K”^JK-i = J. (C7) 

(ii) Let us divide the unitary matrix U into its real 
and imaginary part as 

U = A + iB, (C8) 


where the unitarity enforces the symmetry of AB’’’ and 
the relation AA^ + BB^ = 1. It is straightforward to 
check that the requirement = U0 U* with U G 17(n) 
is equivalent to 



(C9) 


We will now prove the equivalence of this last formu¬ 
lation with M G Sp(2n,R) n 0{2n). If decomposition 
(C9) is assumed, the orthogonality and the symplec- 
ticity of the matrix can be explicitly checked showing 
that indeed MM^ = 1 and MJM^ = J. Viceversa, if 
M G Sp{2n, K) n 0(2n), then 

JM = M'lT J = MJ, (CIO) 


which yields the block structure of (C9). The required 
relations among the blocks follow trivially from the or¬ 
thogonality of the matrix. □ 


Theorem 2. Every symplectic orthogonal matrix W G 
Sp{2n, lR)nO(2n) can be transformed in its standard form 
Wq by means o/M G Sp(2n, R) C 0{2n): 

MTWM = We. (Cll) 


Proof. From Lemma C.l it follows that Wc can be writ¬ 
ten in terms of a unitary matrix U G 17(n), that can be 
diagonalized with a unitary matrix O G U{n). Setting 
Me = O 0 O* one has: 


We 


U0U* = Me 



Mt, (C12) 


with being the eigenvalues of U. Being We = 
K WK“^ by definition, and noticing that 


=K- 


0 , 



0 , 


K, 


(C13) 


setting M = K-i MeK Eq. (C12) yields 


W = MWeMT. (C14) 


The definition of Me and Lemma C.l eventually assure 
that M G Sp{2n, R) C 0(2n). □ 


Appendix D: Minimum in s = 1/2 for two-mode 
squeezing and linear mixing of thermal Gaussian 
states 

In this appendix we want to show that the minimum in 
s, that appears in the definition of the Quantum Chernoff 
Bound (27), is reached for s = 1/2 when the two classes 
of states described in the main text are considered, no 
matter the parameter A used in the optimization set S. 
Similarly to the proof in [33] for the A = 7r/2 case, the 
idea is based on the convexity in s of the Quantum Cher¬ 
noff Bound [21]. Instead of showing the symmetry of the 
involved quantity 

Tr[pV-^], a = UApU\, (Dl) 

under the exchange s GG 1 — s, we will show that its 
derivative is zero on s = 1/2 when the considered state 
is chosen among those obtained from thermal states via 
linear mixing or two-mode squeezing. This will yield a 
statement that holds for every A. 

Using the cyclicity of the trace the derivative of (Dl) 
can be written as 

^Tr [pV^-'*] = Tr [(log p - log a) pV^"'*] , (D2) 

where up to local unitary maps we can always assume p 
as having null first moments. The same hypothesis can 
be made on ct, since the optimization in Ua can always be 
performed among Gaussian operations purely quadratic 
in f (see Sec. Ill for the proof). Being defined as Gibbs 
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state of an Hamiltonian at most quadratic in f, an undis¬ 
placed Gaussian state can always be written as 

p = (D3) 

where Hp is diagonalized into some Dp in terms of the 
same S S Sp{2n) that intervenes in the Williamson de¬ 
composition (14): 

Hp = S-iTDpS-b (D4) 

Therefore also the Gaussian unitary operation Ua, that 
transform the covariance matrix as in (11), alters the 
Hamiltonian in a similar way: 


where p{s), a(s) are the Gaussian states proportional to 
p^, and {•,•}+ is the anti-commutator. In [38] is 
shown that the product of two Gaussian states pip 2 has 
a characteristic function (B9) described by a matrix 

ri2 = -zo -h (r2 -H ifi) (Ti -H r2)”^ (Ti -h zO) , (D7) 

where ri,r 2 are the covariance matrices of pi, p 2 and 
is the symplectic form (5). From definition (7), and 
being 

Tr[pVi-®] = Tr[p5CTi-"p4] e R, (D8) 

it follows that, up to an irrelevant factor, the quantity 
that has to be zero for s = 1/2 is the real part of 


logp-loga = -Hp + H, = ifT f. 

(D5) 

Up to a normalization factor proportional to 
Tr(p^)Tr(cr^“®), the derivative in which we are interested 
in reads: 

^ - Hp)^^ Tr [p(s)a(l - s){r., r ,}+], 

(D6) 


Tr 


(u^'"HpU)(i-Hp) ri2 


(D9) 


where we need to substitute r 2 = and 

Fi = in expression (D7) for ri 2 [as in the main text, 
here F^®) is the covariance matrix of p(s), obtained sub¬ 
stituting its symplectic eigenvalues {vj}j with {As(zzj)}j 
(31)]. 

Exploiting (D4), and labelling as Dg the block diagonal 
matrix defining the thermal part of F^®), the whole real 
part of (D9) can be rearranged as: 


Tr 


D„ 


OTDF^O + D/i, O-iTD/Z^O 


d; 


+ Dp 


1 


1 


OTD,0 + Di_s 0-iTDi_,0 


D, 


, (DIO) 


where the matrix O is defined as S^^U/iS. Notice that 
to prove that such expression is zero when s = 1/2, it is 
enough to show that the quantity 


Tr 



(otDi/2 


1 

O-HD 



(Dll) 


is invariant under the exchange O e-)- 0~^ for every block 
diagonal matrix Dp and Dj ^/2 (whose blocks are multiples 
of the identity). Equivalently, from definition (36) and 
Uy 4 explicit expression (42), this quantity must be an 
even function of A when the symplectic matrix S is chosen 
to be Ssq{r) (46) or SimW (49). This can be checked 
by explicit calculation with the help of a software that 
allows for analytic manipulations, since the high number 
of parameters involved makes it not trivial to obtain by 
hands. 


Appendix E: Relation with Entanglement and total 
number of photons 


GDS for a two-mode squeezed thermal state and other 
quantities: the logarithmic negativity and the total num¬ 
ber of photons. 

Proof of Proposition 1. Eollowing the structure of the 
statement, we divide the proof in two parts: first we 
show that pure states give the minimum GDS allowed 
for a generic state when the entanglement is hxed, and 
then we prove that the whole region above such value is 
actually admissible. 

• We need to prove the first statement, which is: 

>0I?4^)(l,l,£(zzi,z.2,r)/2). (El) 

This can be done in three steps. At first notice 
that the Gaussian DS (59) reaches its minimum for 
a given r when A_ = 0, i.e. when vi = V 2 - In the 
particular case of pure states = 1^2 = !)> this 
yields: 

QVS^^\vi,V2,r) > gVS^^\l,l,r). (E2) 


In this appendix we will prove the two propositions 
stated in sub-Sec. HIB, concerning the relation between 


Secondly, on a pure state GDS is a monotonically 
increasing function with respect to the absolute 
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value of the squeezing parameter r. This can be 
seen deriving Q'DS^^\l, l,r): 


d 

dr 


gvs^^\iA,r) 


sinh(2r) / sin^(A/2) 
[cosh^ r + cot^(A/2)] ^ 


(E3) 


quantity of that has the same sign of r. Therefore 
in order to conclude this part of the proof, we just 


need to show the following inequality: 


2\r\=£{l,l,r)>£{iyi,iy2,r). (E4) 


From definitions (63), (64) and (65) given in the 
main text, an explicit expression for £ can be de¬ 
rived: £ = max{0,/(:/!, ?')}, where 


f{vi,V2,r) = “2 log 


cosh(4r)a:^ -\- x^_ — (cosh(4r)a;^ -|- 


with 


^±=2 ^ ^ 2 ) • 


(E5) 


If f{i'i,i' 2 ,r) < 0 in (E5), the thesis is trivial. 
We need to prove that the result holds even when 
/ > 0. To do so let us suppose without loss of gen¬ 
erality that r > 0. After some manipulations (E4) 
can be shown to be equivalent to the inequality: 

e'*’’ [(21/11/2)^ - {vi + V 2 f] > 

[(*^1 + J^2)^ - 4 ] -b 2(i/i - z^2)^, (E6) 

where both sides are positive. Since the left-hand 
side increases with r, while the right-hand one de¬ 
creases with it, we just need to prove the inequality 
for r = 0, which is the worst-case scenario. In this 
case (E6) can be rearranged as 

I+ ^ 2^2 >^2 ^^ 2 ^ (E7) 

which holds when , 1^2 > I. 

• The second point of the proposition states that two¬ 
mode squeezed thermal states cover the entire area 
above the threshold of pure states and below the 
upper bound = 1, in the QVS^^'^ — £ plot 

(see for example Figure 2). Instead of fixing £, 
varying A on a ^-dependent region, it is easier to 
do the opposite. Moreover, we can consider only 
symmetric states with vi = 1/2 = v- Their GDS 
does not depend on the actual value of v, and it 
is a monotonic increasing function of the squeez¬ 
ing parameter |r|, with values in [0,1[ (see Figure 
1 or Eq. (59) with A_ = 0). The logarithmic neg¬ 
ativity on the other side has a value of 2|r| — log:^ 
when 2|r| > logz^ and zero otherwise. Therefore, 


the parameter 0 < A < 1, defining the value of 
GDS, fixes the squeezing |r|, while £ determines 
the needed v >1, that can be obtained only if £ is 
not too big with respect to the chosen |r| (or equiv¬ 
alently A). It is easy to see that such extreme value 
for £ exactly coincides with the one associated with 
pure two-mode squeezed states. 


□ 

Proof of Proposition 2. As for the previous proposition, 
we divide the proof in two parts: first we prove that pure 
states achieve the maximum GDS for two-mode squeezed 
thermal states when the total number of photon is fixed, 
and then we show that the whole area below such value 
(and above = 0) is admissible. 


• Let us begin by proving 


gVSA^^H^i^^2,r) < 

^l,l,arcsinh . 

(E8) 


Using Eq. (59) for in the considered class 

of states, exploiting relation 


sinh^ 


arcsinh 



N,g{Nsg + 2), 


(E9) 


as well as expression (71) for Ngqfui, 1 ^ 2 , r), one can 
show that (E8) is equivalent to 


cosh^(2r) < cosh^(2r) ^ ^ ^ ^ 




(ElO) 


We are now going to prove the stronger inequality 


obtained from (ElO) by dropping the last squared 
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term on the right, so that the dependence upon 
r disappears. Using definitions (53) for A±, after 
some manipulations the new inequality can be writ¬ 
ten as: 

1 - 1/11/2 < - 1, (Ell) 

which is evidently true being 1^2 > 1- 

• Let us consider the region of the plot — Ngq 

delimited by the boundary of the pure states and 
the axis — 0. We need to prove that for 

every couple {Nsq,A) in such region, there exists 
a two-mode squeezed thermal state with number 
of photons Nsq and GDS equal to A. To do so we 


can consider once again only symmetric states with 
vi = 1/2 = v. Moreover, instead of fixing Ngq while 
moving A on a As,j-dependent region, we can do 
the opposite. Indeed, for a symmetric state GDS 
is a monotonic function of |r|, which is therefore 
fixed by A. The constraint given by imposing a 
total number of photons Ngq then reads: 

i/cosh(2r) - 1 = (E12) 

that can be solved for v > 1 only if Ngq is big 
enough. It is easy to see that such extreme value 
for Nsq exactly coincides with the one associated 
with pure two-mode squeezed states. 

□ 




